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ABSTRACT 


A  brief  review  is  given  of  attempts  to  formulate 
a  unitary  field  theory  involving  the  electromagnetic 
field.  This  thesis  concerns  itself  with  another  such 
attempt,  based  on  a  proposal  by  Dr.  H.  Schiff  for  non¬ 
linear  six-vector  field  equations. 

The  static  case  of  these  equations  has  been  con¬ 
sidered  and  analytic  solutions  were  obtained  for  two 
special  forms  of  the  static  equation.  A  third  form  of 
the  static  equation  was  solved  on  the  LFP-30  computer 
and  studied  in  some  detail.  It  is  shown  that  the  latter 
solutions  describe  particles  with  no  charge. 

For  the  third  form  of  the  equation,  three  discrete 
solutions  were  found  on  the  computer.  Using  a  heuristic 
expression  for  the  total  energy  of  these  particles,  a 
value  was  obtained  for  the  coupling  constant  of  this 
field  which  was  of  the  order  of  the  coupling  constant  for 
"strong"  interactions. 

Interaction  potentials  between  the  particles  were 
evaluated,  using  analytic  approximations  to  the  computer 
solutions.  A  ratio  of  the  energy  of  the  first  particle 
to  that  of  a  compound  particle  composed  of  the  first  two 
particles  was  found  to  be  comparable  to  the  ratio  of 
pi-meson  to  nucleon  mass. 
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1. 


I.  INTRODUCTION 

Maxwell* s  classical  electromagnetic  equations  are 
dualistic  in  form.  To  describe  the  physical  world  two 
separate  entities  are  needed:  fields  and  particles.  The 
particles  form  the  sources  for  the  field  and  interact 
with  the  field  but  are  not  a  part  of  the  field.  They  are 
assumed  to  be  point  particles  with  their  only  character¬ 
istics  being  charge  and  mass.  The  dualistic  nature  of 
the  theory  is  clearly  exhibited  by  the  Lagrangian  for  a 
particle  field  system:1 

S.  -  f fJ  d  t  (1.1) 


in  which  the  field,  the  particles  and  the  interaction 
are  introduced  separately.* 

The  electromagnetic  field  tensor  is  defined 
as  the  curl  of  a  four-vector  potential,  Aj_,  in  the 
following  manner: 


Fiw  =  3Ak  _ 


(1.2) 


Because  in  the  classical  theory  the  four-potential  is 
introduced  essentially  as  a  mathematical  artifice,  one 


See,  for  example,  Landau  and  Lifschitz,  THE  CLASSICAL 
THEORY  OF  FIELDS,  ( Addison-Wesley  Press,  Inc.)  p.  77. 

*  See  Appendix  for  definitions  of  symbols. 
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of  the  requirements  which  has  been  imposed  on  the 
classical  theory  is  that  of  gauge  invariance,  i.e.  the 
equations  are  to  be  invariant  under  the  transformation, 


A;  -  Ac  - 


is 

9*1 


(1.3) 


where  S  is  some  arbitrary  scalar  function. 

The  total  momentum  and  energy  of  the  electro¬ 
magnetic  field  can  be  expressed  as  the  integrals  of 

a  tensor  T.^  called  the  energy-momentum  tensor.  These 

2 

expressions  are: 


p„ 


»,  3  = 


1  ^  ) 


(1.4) 


The  integration  is  extended  over  a  hypersurface  containing 
all  of  three  dimensional  space.  If  the  integration  is 
carried  out  on  the  hyperplane  t  -  a  constant,  we  then  have: 

P r-'  ~  ”  T J  V 
£  =  -  C TA4d v 

(1.5) 

2  Landau  and  Lifschitz,  loc .  cit.,  P.  82-86. 
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3. 

Because  of  the  dualistic  nature  of  the  theory, 
profound  unsolved  problems  exist  which  emphasize  the 
need  for  a  more  satisfactory  theory.  Some  of  these 
difficulties  are  considered  below. 

In  the  spherically  symmetric  static  case,  the 
electric  field  for  the  electron  follows  the  well-known 
Coulomb’s  law, 


Using  eqns.  (1.4),  (1.5)  and  (1.6),  we  obtain  the 
following  expression  for  the  total  energy  of  the 
electron . 


£ 


f*‘>v 


(1.7) 


If  the  electron  is  a  point  particle,  the  integral  must 
be  taken  up  to  the  origin  and  the  total  energy  becomes 
infinite . 

Lorentz  and  Abraham  removed  this  difficulty  by 
assuming  a  finite  size  for  the  electron.  If  the  electron 
is  assumed  to  be  a  sphere  of  radius  "a”,  the  integral 
can  be  evaluated  and  the  total  energy  becomes, 

_e_2~ 

2  CL 

(1.8) 

This  theory  introduces  the  necessity  for  some 
cohesive  force  to  hold  the  electron  together,  which 
must  be  nonelectromagnet ic  in  origin. 
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4. 

Another  difficulty  with  the  finite  electron  in 
the  classical  theory  is  that  the  four-momentum  does 
not  obey  the  relativistic  transformation  laws  for  a  four- 
vector.  We  consider  an  electron  moving  with  a  velicity 
"v"  along  the  z-axis .  Then  using  eqn.  (1.5),  v/e  have 
that 

£  =  -  f  T44dV 


(1.9) 


in  the  rest  frame  of  the  electron.  In  the  frame  of  the 
observer,  we  have 

£'«  -  [T44'  d\l’ 

Pi  -  [T3/dV' 

(1.10) 


Since  in  the  rest  frame  of  the  electron  T-^  is  zero, 

the  result  of  using  the  relativistic  transformation 

3 

formulas  is 


T 


3  3 


dv1^  y  1- y.*  av 

V  •  c  *• 


(1.11) 


3  Landau  and  Lifschitz,  loc .  cit.,  p.  18. 
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5. 


Substitution  of  the  above  expressions  into  eqn.  (1.10) 
gives: 


-*.*3  jV,,dv 


(1.12) 


Since  is  a  four-vector,  a  Lorentz  transformation 
must  give: 


£ 


Ta4  clV 


V 

C  'Z- 


r 


A 


_  v3 
c  *■ 


—  )  4  4-  ^  \/ 


(1.13) 


If  eqns.  (1.12)  and  (1.13)  are  to  be  equivalent,  the 
following  equation  must  be  satisfied. 


f  T„ ,  d  V  =  o 


(1.14) 


By  a  similar  argument,  one  can  show  that  the  integrals 
of  T-q  and  T22  in  the  rest  frame  of  the  electron  must 
also  be  equal  to  zero* 


6 . 

For  an  electron  with  a  finite  radius  "a",  these 
integrals  are  not  equal  to  zero.  Their  values  are: 

[TuJV  -  f T” av  -  -  e~  %: 

(1.15) 

Thus  the  finite  electron  does  not  have  the  proper  trans¬ 
formation  properties. 

Aside  from  the  "dualistic"  difficulties,  one  also 
recognizes  limitations  of  the  classical  theory  vis-a-vis 
quantum  processes.  One  example  of  this  is  the  phenomenon 
of  photon-photon  scattering.  Since  Maxwell's  equations 
are  linear,  the  principle  of  superposition  holds.  Thus 
in  the  classical  theory,  one  electromagnetic  wave  could 
not  interact  with  another  one  and  be  scattered  by  it. 

Considerations  of  this  sort  have  led  to  attempts 
to  abandon  the  dualistic  viewpoint  and  adopt  a  Unitarian 
one,  in  which  the  electromagnetic  field  and  matter  are 
considered  to  be  different  states  of  the  same  physical 
entity. 

In  this  thesis  we  will  examine  attempts  to  formu¬ 
late  a  unitary  theory  in  which  matter  is  considered  as 
essentially  electromagnetic  in  nature.  This  leads 
naturally  to  a  non-linear  generalization  of  Maxwell's 
equations . 
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7. 


In  Section  II  the  theory  of  Mie  is  outlined  briefly 
and  its  limitations  are  indicated.  In  Section  III  we  con¬ 
sider  the  theory  of  Born  which  attempts  to  overcome  some 
of  the  difficulties  inherent  in  Mie’s  theory  but  which 
has  unsatisfactory  features  of  its  own.  In  Section  IV 
we  give  an  example  of  another  attempt  by  Rosen,  which 
requires  the  introduction  of  a  subsidiary  field.  In 
Section  V  we  examine  in  some  detail  suggestions  for  a 
unitary  theory  proposed  by  Dr.  H.  Fchiff  in  which  the 
quantum  mechanical  field  equations  provide  a  point  of 
departure  for  the  construction  of  a  non-linear  six-vector 
theory.  Finally  in  Section  VI  certain  tentative  con¬ 
clusions  are  drawn  with  regard  to  this  latter  theory. 


8. 


II.  MIE'S  THEORY 


The  first  attempts  to  obtain  a  non-linear 

generalization  for  the  electromagnetic  field  equations 

L. 

were  made  by  Gustav  Mie  .  Mie  assumed  that  the  electro¬ 
magnetic  field  is  the  fundamental  entity  and  that 
matter  is  just  a  manifestation  of  this  field.  In  this 
picture  the  electron  is  not  a  foreign  body  in  the 
electromagnetic  field  but  a  special  form  of  the  field. 

It  is  not  a  point  charge  but  a  distributed  charge  with 
most  of  the  charge  concentrated  at  the  center.  The 
principle  of  relativity  was  also  assumed  to  be  valid. 

Mie  describes  the  electromagnetic  field  entirely 
in  terms  of  the  so-called  entities  of  quantity;  D,  H, 
p  ,  and  7;  and  the  entities  of  intensity;  E,  H,  and  7f. 
These  entities  obey  the  equations: 


Curl  H 


_J_  3  V  c  A  TT  J 

c  2-t  C 


d  iv  V  =  4lr  P 


Curl  A 


(2.1) 

(2.2) 

(2.3) 

(2.4) 


4 


Annalen  der  Physik,  Vol  32?  P*  5H  (1912);  Vol.  39. 

p.  1,  (1912). 
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9. 

These  equations  can  be  put  into  four-dimensional 
notation.  The  entities,  Fik,  P-k,  and  Aj_  are  defined 


(2.5) 


■*1,1,3  -  *  >  H  1  *  J 

X  4  - 

let 

(2.7) 

J4  = 

1  C  p 

(2.8) 

A  Ij  2,  3  =  Ax,  yj  4  3 

A* 

(2.9) 

The  Latin  indices  take  on  the  values  from  one  to  four, 
and  all  repeated  indices  are  summed  from  one  to  four. 
With  this  notation  eqns.  (2.1)  -  (2.4)  can  be  written, 

3  -  4  TT  j  c 

3  x  k  c 


(2.10) 


. 


. 

(  „  )  • 

f  •  j  ro '  ••• '  •  * 

.  - :  •  o '  o 

,  (  ' r  i  -  r  „  . 


10. 


r  u  ^  (2.11) 

The  entities  E  and  B"  are  the  electric  and  magnetic 
intensities,  respectively,  and  "D  and  H  are  the  electric 
and  magnetic  field  strengths.  The  entities  "A  and  <-P  are 
the  vector  and  scalar  potentials,  and  J  and  p  are  the 
current  and  charge  densities. 

Two  more  equations  follow  from  eqn.  (2.10)  and 
(2.11) .  These  are 


-  O  (2.12) 

”3  x  i 


ifU  +  3fLl  -t-  9  Fj. —  o 

D  x  £ 


(2.13) 


These  equations  can  be  written  in  three  dimensional 
notation  as: 


d  \  v  J  +  ^  9 

3-t 


(2.12a) 


(2.13a) 


curl  ■+■  -£■ 


9  B  _ 

9t 
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(2.13b) 
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11. 


The  entities  of  quantity,  P.k  and  j  are  complex 


functions  of  the  entities  of  intensity,  Fik  and  A 


This  relation  is  found  through  a  variational  principle. 
Mie  considers  a  "world  function"  or  Lagrangian  which  is 
a  function  of  the  entities  of  intensity.  Then  the 
equations  of  motion  are  such  as  to  make  the  space-time 
integral  of  this  Lagrangian  a  minimum. 


S  j"  jC  C  F-w,  Ai)Jt=  o 


-  r  ^  5/1  '  n  f  ^  =  o 


If  the  fields  go  to  zero  at  large  distances,  the 
divergence  theorem  can  be  used  to  transform  this  integral 


to  the  form: 


<1 


. 

1 


- 


12. 


Since  the  variation  of  the  A.^  is  arbitrary,  the  integrand 
must  be  identically  zero.  Thus  we  have 


3  3  <=£  4  3  c£  —  q 


(2.14) 


If  this  equation  is  to  correspond  to  eqn.  (2.10),  the 
following  relationship  between  the  entities  of  quantity 
and  the  entities  of  intensity  must  be  satisfied: 

P.k  =  -  ^_sL  (2.15) 

3  Fck 


(2.16) 


Of  course  the  relationship  is  maintained  if  the  right 
hand  sides  of  the  eqns.  (2.15)  and  (2.16)  are  multiplied 
by  an  arbitrary  constant.  The  constant  is  chosen  so  that 
for  Maxwell's  Lagrangian,  <£  =  -  -fab  f\  £  ,  we  have 

Piw  -  Fck 


In  three  dimensional  notation  these  equations 


become : 

hf  =  —  4  TT 

3  e 

"Q  -  4  ti  S 

3  E 


(2.15a) 


3  <£ 

-  9  ^ 

3  </> 


(2.16a) 


where 


13. 


3  /  = 

9-1? 


a 


3  of 

-5  0, 


*  J 


■3Bk 


l 


The  entities  of  intensity  can  also  be  considered 
as  complex  functions  of  the  entities  of  quantity.  This 
relationship  is  determined  through  a  Hamiltonian  which 
is  a  function  of  the  entities  of  quantity.  The  Hamiltonian, 
w-  ,  is  defined  by  Mie  as: 


5/  =  __L  Fa  Pa 


9TT 


Ac  jt  -  o£ 


(2.17) 


Using  eqns.  (2.15)  and  (2.16),  we  have, 

d  ~  sv  Fik  d  PiK  ^  4: 

Using  the  chain  rule  for  differentials,  we  obtain  the 
following  equations. 


24  -  -  f,, 

3Rk  4  ti 

3—^  -  A_h 
c 


(2.18) 


(2.19) 
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In  three  dimensional  notation  these  equations  become: 

71  ^  ^ 


B  -  -  4 


3T7 


r  -  4  ti  ^ 


(2.18a) 


l- 


^  w_ 
3? 


c_p  -  —  3  ^ 


(2.19a) 


By  using  a  variational  principle  on  the  Lagrangian, 
we  obtain  eqn.  (2.10)  from  eqn.  (2.11).  We  can  reverse 
this  procedure  by  carrying  out  a  variation  of  the  Hamiltonian, 
starting  with  eqn.  (2.10),  and  using  eqns.  (2.18)  and  (2.19) 
as  definitions. 

The  variational  principle  is  : 

s  y  xjT=o 


3  £  Pc ^  j  <1  )  clx  :  0 

^Rk 


Using  eqn.  (2.10),  we  have 


S  j  C  -  Sz_  _3__  ^  ^  W 

4r  3Xk 


[  Fwlt  SP.w  -  )  d-t  =  o 


15. 


The  divergence  theorem  gives: 


-  -jr  Rk  ^  P«.k  — 


5 SR.w')d'c=o 


iJi  s  Pew  = 


I.  (  5-Ai 

2  v  3xk 


=£*  'i  s  Rk 
) 


Therefore , 


=  G 


and  F.  k  -  -  — - i 

Thus  minimizing  the  space-time  integral  of  the  Hamil¬ 
tonian  leads  to  eqn.  (2.11) . 

The  problem  of  formulating  a  unitary  field  theory 
has  thus  been  reduced  to  one  of  determining  the  correct 
form  of  the  Lagrangian  or  Hamiltonian.  Although  their 
actual  form  is  not  known,  the  principle  of  relativity 
requires  that  they  be  four  dimensional  scalars  or  invar¬ 
iants.  Therefore,  they  must  be  functions  of  the  invariants 
which  can  be  constructed  from  the  field  entities. 

It  can  be  shown  that  there  are  only  four  indepen¬ 
dent  invariants  which  can  be  constructed  from  the  entities 
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16. 


of  intensity,  to  be  used  in  the  Lagrangian.  These  four 

are*  F  ik’  FikF*ik’  ai  ’  and  ^ik^k^2*  The  dual  tensor, 
F*ik?  is  defined  as  ieikimF^m,  where  e^lm  is  the  com¬ 
pletely  antisymmetric  pseudo-tensor  of  the  fourth  rank. 


eiklm 


iklm  an  even  permutation  of  123^ 
iklm  an  odd  permutation  of  123^  (2.20) 

any  two  indices  the  same 


Therefore 


(2.21) 


The  invariants  in  three  dimensional  notation  are: 

«  2cb*-e>) 

Pih  Ft*  =  -M  L  E) 

A.1  «  A2-  (t* 

(RwA,)2=  Lj 


(2.22) 


17. 


Similarly,  there  are  only  four  independent  invar¬ 
iants  which  can  be  constructed  from  the  entities  of 
quantity.  These  are: 

PU*  .  2  c h  d *) 

PlU  P.l,**  -4a  Cl7  73) 

j‘  =  J‘-CV  (2.23) 

CP.UO 2  -  (?**  *cfl5)1-  (Td)z 

A  conservation  law  can  be  drrived  using  the  field 
equations  (2.10)  and  (2.13).  If  the  Lagrangian  is  differen¬ 
tiated  with  respect  to  the  coordinates,  the  result  is: 

:2j£  ..  _L  2j£  3  u  -+  2J  rLii 

3  *  t-  ^  3Pi,j g 

Using  eqns.  (2.15)  and  (2.16)  we  have, 

^  cT  _  __  _ i _  ^  /  3  F~w  i  \  +■  J  i_  ^ 

3xt'  £71  ^  3xi  i  c  9  *  i 

Uwing  eqn.  (2.13), 


3  dC 

1 

^WJ.  ( 

+  3Fik 

)  +  A* 

3  At 

d  X  £ 

*3  T7 

^  3  x  i 

■L 

3  *XjL 

3^-c 

L 

Ph* 

3F|; 

-+ 

J* 

3  /4jt 

4  77 

c 

3  X(,' 

-  1 

4  T7 

9X  ^ 

CPki 

Sc) 
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Using  eqns.  (2.10)  and  (2011), 

dS  ^ 


-i-  (FcxRu  )  +-/AA 

4TT  3Xk  l  ?Xi 


9* 


^  )H  -t  i±  2Aj 

«-  /  C. 


”  ~r  -2-  t  Rj.  Ph  j  )  +  jk  3/1 ; 

471  3i>v  -s:  3*,, 


Using  eqn.  (2.12) , 

-i,  SiK<C  ,  - 


3 1. 


-J—  3_  tF;i  Pa)  +  _3_  Acj* 

4-  "n  9  h.  3xk 


where  is  the  Kronecker  delta.  If  the  tensor  T. 


ik 


is  defined  as: 


"Tew  -  $:k  £  +  — i —  F;i  Rzi  -  4 


4tt 


c  J  k 


(2.24) 


Then 


^  ~F  i  k  =  o 
9Xh 


Therefore  the  integral  of  this  tensor  over  a  hypersurface 
containing  all  of  three  dimensional  space  is  conserved. 
This  tensor  is  called  the  energy-momentum  tensor.  The 
total  momentum  or  energy  of  the  field  can  be  found  by 
using  eqn.  (1.4)  . 


, 

0 

;  " 

t 

<* 

0 

■  •> 

.  ) 


n  ■  '  '  ‘  '  '  '  ;  i? 

"i  ' 


•  ■  \ 

iijr 

19. 

Although  Mie  has  developed  a  general  form  for 


the  electromagnetic  field  equations  through  the  use  of 
a  Lagrangian,  he  gives  no  way  of  determining  the  form 
of  this  Lagrangian.  His  only  criteria)!  are  that  it 
be  a  relativistic  invariant  and  that  at  large  distances 
from  the  center  of  the  electron  his  Lagrangian  should 
approach  Maxwell's  Lagrangian. 

Mie  states  that  since  the  results  of  applying 
his  theory  would  be  compared  with  experiments  for  the 
static  case,  he  uses  only  the  invariants  Fj^  and  A ^ ? 

since  in  the  static  case  F^f*^  is  equal  to  zero  and 
2 

(FifcA^)  can  be  expressed  as  the  product  of  the  first  two 
invariants.  He  has  considered  a  Langrangian  of  the  form: 


Then, 


He  has  given  the  solutions  for  n  -  3*  He  has  shown  that 
in  this  case,  one  obtains  a  highly  concentrated  charge 
distribution,  but  it  is  unstable  in  the  field  of  another 
electron. 
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2 

The  use  of  the  invariant  A ^  in  Mie*s  theory, 
introduces  the  potentials  directly  into  the  field  equations. 
Therefore  the  equations  of  motion  are  not  gauge  invariant, 
i.e.  his  theory  depends  on  the  absolute  value  of  the 
potentials.  Thus  in  the  static  case,  if  the  external 
potential  were  changed  by  a  constant,  the  form  of  the 
solution  could  change.  Also  since  the  equations  are 
derived  on  a  classical  basis,  there  is  no  possibility  of 
explaining  the  spin  of  the  electron. 

Despite  its  shortcomings,  Mie*s  general  formalism 
has  been  the  starting  point  for  most  unitary  field  theories, 
such  as  the  one  by  Born,  which  we  consider  in  the  next 
section. 
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III.  BORN'S  THEORY  OF  THE  ELECTRON 

In  193^i  Born^}  j_n  colloboration  v/ith  Inf  eld, 
developed  a  Lagrangian  to  describe  the  electron.  Their 
theory  was  gauge  invariant  in  that  only  the  fields  and 
not  the  potentials  appear  explicitly  in  the  Lagrangian. 
They  developed  a  Lagrangian  such  that  the  action  integral, 
j'  d 't  ,  was  an  invariant  under  the  transformations 
of  general  relativity. 

Using  a  covariant  tensor  a^-p  the  simplest  form  of 
the  Lagrangian  so  that  the  action  integral  is  invariant 
is  the  determinant  J |akl|  .  The  tensor  akl  is  split  into 
a  symmetrical  and  an  antisymmetrical  part. 

The  symmetrical  part  is  identified  with  the  metric  tensor 
and  the  antisymmetrical  part  with  the  electromagnetic 
field  tensor. 

The  simplest  Lagrangian  which  can  be  constructed 
with  these  entities  is: 

+  /"  I  *3  h  1 1  +  *  /lK7 

5  Proc.  Roy.  Soc.  London,  A 144,  p.  425,  1934. 
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The  last  term  can  be  neglected  since  it  gives  no 
contribution  to  the  field  equations.  In  cartesian  coordin¬ 
ates  and  with  the  notation  of  special  relativity,  we  have, 

<£  '  K  (  /  I  +F-G1-  '  )  (3.1) 

F «  Rk  =  Bx-F'  (3.2) 

TV-  b1 

G  (3.3) 

The  constants  are  chosen  so  that  in  the  limit  of  weak 
fields,  this  Lagrangian  reduces  to  Maxwell’s  Lagrangian. 
Here  "b"  is  a  constant  with  the  dimensions  of  field 
strength. 

This  Lagrangian  sets  an  upper  limit  "bM  on  the 
strength  of  the  fields,  just  as  the  transformations  of 
special  relativity  set  an  upper  limit  "c"  on  the  velocity 
of  a  particle. 

The  entities  of  quantity  can  be  found  by  using 
eqns.  (2.15)  and  (2.16).  Thus  we  have, 


Few  -  i  G  Rk 

(3.4) 

J I  ■+  F  -  G 

a  ° 

(3.5) 

H  -  B*  -  -fr  (B-g)  E 

J  l  +  F  -  G 2 


(3.6) 
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D  -  E  -  ~b  (  B  •  E  )  B 
/  /  +  r-  GL 

Then  the  equations  of  motion  become: 

3^  ,o 


(3.7) 


(3-8) 


Cur 


H  -  -L  £D 

C  9t 


-  O 


(3.9) 


i  v 


D 


(3.io) 


By  using  eqn.  (2.24),  the  energy-momentum  tensor, 
Tj^,  can  be  determined.  For  this  Lagrangian,  the  energy- 
momentum  tensor  is: 


/l  +  F  -  6  2 


-  by  ,/m  F-C  2  I 

(3.11) 


Born  and  Inf eld  have  considered  the  static  case 
with  central  symmetry.  In  this  case  the  fields  become: 


d  i  v  D  -  o 

D  r  =  e. 

r2- 


(3.12) 


e_ 
r  "*■ 
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E 


r 


bv4 


(3.13) 


qc> 


x/I+J* 


ro>/-|  (3.14) 


The  potential  is  finite  and  continuous  everywhere  and 
goes  as  e/r  asymptotically.  The  value  of  the  potential 
at  the  origin  is  1.854e/r0.  The  electric  field  is  also 
finite  everywhere  but  it  has  a  discontinuity  in  direction 
at  the  origin. 

Born  and  Infeld  consider  the  divergence  of  "E  as 
defining  a  spatial  charge  distribution,  which  they  call 
the  "free  charge"  density.  The  actual  charge  density 
is  zero  since  the  potentials  do  not  appear  explicitly 
in  the  Lagrangian.  For  the  static  case,  this  free  charge 
density  is: 


_ _ e _ _ 

C=  jnr.«r  [.+  (-fJ4J3/i 


(3.15) 


It  is  evident  that  the  charge  density  becomes  very  small 
0.  The  total  free  charge  over  all  of  space  is: 


e 


when  r  >  r 
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The  value  of  the  constant  "b"  in  these  equations 
is  determined  by  calculating  the  total  energy  of  this 
particle  and  equating  this  to  the  rest  energy  of  the 
electron.  For  the  static  case  we  have, 


(3.16) 


Using  eqn.  (1.5),  we  have, 


If  eqn.  (3.13)  Is  substituted  into  the  above,  the  result 
is: 


£  =  1.2361  e_!  =.  rr  C  1 

Co 


(3.17) 


where  "m"  is  the  rest  mass  of  the  electron.  This  gives 

To  ~  2.2gX/o-13  cm  , 


b  =  9.16  x  / 6lS  e.s.v. 


(3.18) 


Because  of  the  small  value  of  the  characteristic  length 
of  the  electron  and  the  large  value  of  the  maximum  field 
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strength  in  this  theory,  experimental  deviations  from  the 
classical  laws  would  not  be  readily  apparent. 

Although  Born's  theory  has  the  desirable  property 
of  gauge  invariance,  it  is  not  truly  unitary,  since  the 

6  *7 

definition  of  charge  is  artificial.  Frankel  and  Feenberg' 
have  shown  that  the  field  equations  must  be  supplemented 
with  equations  of  motion  for  the  electron  for  a  complete 
description. 


6  Proc.  Roy.  Soc.  London,  A146,  p.  933?  193^+  • 
Phys.  Rev.,  Vol.  47,  p.  148,  1935. 
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IV.  ROSEN'S  FIELD  THEORY 


In  this  section,  we  indicate  briefly,  a  somewhat 

o 

different  formulation  due  to  Rosen.  He  has  shown 
that  any  Lagrangian  which  contains  only  the  fields  ex¬ 
plicitly,  has  either  a  singular  point  or  a  discontinuity 
in  direction  of  the  field  at  the  origin  in  the  spherically 
symmetric  static  case.  This  requires  that  the  field 
equations  be  supplemented  with  equations  of  motion  for 
the  particles.  For  a  truly  unitary  theory,  the  potentials 
must  appear  directly  in  the  Lagrangian.  However,  if  this 
were  the  case,  the  equations  would  not  be  gauge  invariant. 

To  remove  this  difficulty,  Rosen  introduced  an 
additional  scalar  function  which  changes  under  a  gauge 
transformation  in  such  a  manner  that  the  total  Lagrangian 
is  gauge  invariant.  He  chose  a  "matter"  function,  Lp  , 
which  transforms  in  the  following  manner: 

L|/ '  -  e  *  £  5 

as  (4.1) 

A  k  Ak  +  ii 

3  *  k. 

Here,  £  is  a  constant  with  the  dimensions  of  inverse 
charge.  Two  simple  functions  of  the  potentials  and  this 
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Fhys .  Rev.,  Vol.  55?  p.  94,  1939. 
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scalar  function  can  be  constructed,  which  are  both  Lorentz 
and  gauge  invariant.  These  are: 


^  4'* 


(4.2) 


■Xk-K 


(4.3) 


^  _  x  e  /4k  QJ  (4.4) 

where  (*)  stands  for  complex  conjugate. 

Rosen  has  chosen  a  Lagrangian  of  the  form: 

<£  =  -  Sii  -  "Xt  +  °"z  ^  (4.5) 

/6TT 

where  <r  is  a  constant  with  dimensions  of  inverse 

length.  Interestingly  enough,  about  ten  years  previously, 

9 

the  same  Lagrangian  had  actually  been  propsed  by  Gordon  , 
Klein10,  and  Schr oedinger11’ 12 ,  although  Rosen  was 
apparently  unaware  of  this. 

9  Z.  Phys.,  40,  117,  (1926) 

10  Z.  Phys.,  32,  895,  (1926) 

11  Ann.  der  Phys.,  (4),  8l  (1926);  82  (1927) 

12  Nature,  Vol.  169,  P.  538,  (1952) 

^8  Phys.  Rev.,  Vol.  62,  p.  436,  (1942) 
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In  using  the  variational  principle,  the  potentials 
and  the  matter  function  are  varied  independently.  This 
gives  an  additional  set  of  equations  besides  the  usual 
Maxwell’s  equations.  For  Rosen's  Lagrangian  the  equations 
are : 

( ~  ^h)  U/  -f  a~  z  ^  =  o 

l  /  ‘  (4.6) 

3  R 1  -  4 tt  j  n  (4.7) 

_L  j*.  =  -1.  4e  (qj*Xk  -  ^  Xt  )  (4.8) 


1*5 

Menius  and  Rosen  J  have  solved  the  spherically 
symmetric  static  case  numerically,  and  have  been  able  to 
obtain  electron-like  solutions.  The  solutions  however 
had  the  undesirable  property  of  having  the  total  energy 
negative . 
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Phys.  Rev.,  Vol.  62,  p.  436,  (1942). 
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V.  FOUNDATIONS  OF  A  UNITARY  FIELD  THEORY 

In  this  section  we  consider  a  proposal  of  Dr. 

H.  Schiff  for  a  unitary  six-vector  field  theory. 

Since  the  best  description  yet  obtained  for  a 
particle  is  provided  by  the  Dirac  equation  for  the  electron, 
it  seems  natural  to  use  this  equation  as  a  basis  for  the 

development  of  such  a  theory,  Dirac’s  equations  can  be 

...  14 

written 


JL.  3  (fi  —  m  4-  CT*  •  p  X 

1 Tt 

(5.D 

1.  ax  .  -  m  X  +  O 4> 

(5.2) 

with  ft  -  c  =  1 . 

The  functions  and  X  are  first  rank  spinors  which 
under  an  inversion  of  the  coordinate  system  transform  as: 


<J>'=  4 

")(/  =  -X  (5-3) 
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Akhiezer  and  Berestetsky,  Quantum  Electrodynamics, 
(Consultants  Bureau,  Inc.)  p.  55. 
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In  eqns.  (5.D  and  (5.2),  p  is  the  momentum  operator, 
-i  grad,  and  cr  is  the  set  of  Pauli  matrices, 


(5.4) 

If  eqns.  (5.D  and  (5.2)  are  added  and  substracted 
the  result  is: 

*._s2_C</>4*)s  m(4-%)  (4  +~A) 

3 1’ 

(5.5) 

=  m  £<£  +"X  -  6) 

Under  spatial  inversion,  (  <b  +  TO  and  (ck-7')  transform 
into  each  other . 

We  now  consider  the  complex  field  vector  •*  H  +•  iE*, 
where  1?  is  an  axial  vector  and  E  is  a  polar  vector.  It 
is  known  from  the  work  of  Laporte  and  Uhlenbeck1^  that 
appropriate  linear  combinations  of  the  components  of 
F  transform  as  the  components  of  a  rank  two  spinor.  Fur¬ 
ther  by  contraction,  one  can  always  form  a  first  rank 
spinor  by  combining  a  second  rank  spinor  with  a  first 
rank  spinor.  One  can  then  show  that  if  B  is  an  arbitrary  first 
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Laporte  and  Uhlenbeck,  Phys .  Rev.,  Vol.  37,  p.  I38O  (193D* 
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rank  spinor,  then  (d*-F  )B  is  also  a  first  rank  spinor. 
That  )B  is  a  first  rank  spinor  will  be  evident  when 

we  consider  the  two  component  form  of  the  Dirac  equation. 

We  now  note  that  under  inversion  of  corrdinates, 

(  6s  •  F)  and  (  o*  •  F*)  transform  into  each  other.  If  we 
now  choose  two  first  rank  spinors  B  and  C  which  transform 
under  inversion  as  ^  and  c £  -  X  respectively, 

then  it  follows  that  (  a*  *  "F)B  and  (  *  F*)C  also 

transform  as  <£  +  7C  and  respectively. 

Eqns .  (5*5)  and  (5.6)  can  be  easily  combined  to 
obtain  a  second  order  two  component  equation, 

=  °  <!•?> 

and  a  similar  equation  for  -  r  .  In  the  presence 
of  an  electromagnetic  field,  this  equation  becomes: 

f(x-sL  -  e?'?jc(D+^)=o  (5.8) 

It  is  now  apparent  that  since  eqn.  (5.8)  is  covariant 
that  (  <?  •  ~F)  (  +  "X  )  must  transform  as  4>  +~X  . 

We  now  replace  the  spinor  (j>  +X  in  eqn.  (5.8) 
by  (  Q*  •  1f)B,  with  the  s  in  the  equation  being  the 

— 9 

potentials  for  F„  Then  the  equation  becomes: 

[(a-3-  B=o  (5.9) 
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This  is  a  relat ivistically  invariant  equation  for  a  self 
interacting  field  of  spin  The  constant  "g"  in  the 
equation  has  the  units  of  charge  and  gives  a  measure  of 
the  interaction. 

Eqn.  (5.9)  suggests  that  for  the  boson  case,  we 
consider  the  following  equation: 

(5.io) 

where  the  operator,  (  L  ~^xn~{r  $  ^  1  ^  ^  is  simply 

the  Klein-Gordon  operator  with  a  self  interaction. 

If  the  constants,  and  c,  are  restored,  the 
equation  becomes 

[  (jL*^  ■+  _S_  m2cJ  1  F  =  0 

L  ^  G  J  (5.1D 

Equation  (5.11)  is  not  gauge  invariant  since  the 
transformation 

Aw'-  Aw  + 

F  '  -  F  e  ‘ 


gives 


(5.12) 


; 


If  we  assume  that  H  and  E  are  real  fields,  eqn. 
(5.H)  can  be  split  into  real  and  imaginary  parts.  Then 
the  field  equation  can  be  written  as: 


-_S?  Ak1- 


rv\ 


+  JL±  A* 

c.  <S>Zh. 


+  i  >U  Fj.  r*\  K  ■=.  O 


(5.13) 


where  F^m  and  F*lm  are  defined  in  eqns.  (2.5)  and  (2.21). 

The  static  case  of  this  equation  is  considered  in  the  next 
section. 

Attempts  to  find  a  Lagrangian  for  this  theory  have 
so  far  been  unsuccessful.  Presumably  in  order  to  obtain 
a  Lagrangian  which  goes  over  into  the  Maxwell  Lagrangian 
as  "m"  and  "g"  go  to  zero,  one  would  require  a  first  order 
boson  equation.  However  since  there  are  a  number  of 
different  types  of  first  order  boson  equations,  all  of  which 
would  lead  to  a  Klein-Gordon  type  of  second  order  equation, 
and  also  since  we  will  restrict  ourselves  to  the  static  case 
we  will  use  eqn.  (5.11)  and  a  heuristic  expression  for 
the  energy. 


I 


rf 


) 


<  i 


■) 


1 


-)  1  I 


1  'I 


L  ,  - 


o 


[  ,r 


i  ’  f>  '  S  i¬ 


ll 


v  ■  ■  -  - 


‘ ;  o 


•ofvtc  ' 


■  ^  •  v  '  i  ••  v  •  ■  ■ : 

' :  "  .  ■  v '  n  •  - 


~0 


■  ■  • 


r  * 


Ci’ 


o 


35. 


VI.  ANALYSIS  OF  THE  STATIC  CASE 
For  the  static  case  and  with  7T  =  0,  eqn.  (5.13) 

becomes 

V2E-  rr\xcx  £1  c^lT-o 

.1**0 

— ^ 

Putting  E  -  -grad  W3  we  get, 

V  2  V  Lfi  -  V  +  -4-1  <-P  z  VC#  -  o 

Jl  l  Jc.  it1 


i  .e 


from  which  it  follows  that  the  expression  in  brackets  is 
a  constant  (Y,  say),  or, 


C7X  (_p  ^  <Z)  -  44  -f  Y' 

^  3 


(6.1) 


It  was  not  possible  in  general  to  obtain  analytic 


solutions  to  this  equation  and  numerical  methods  were 
necessary.  However  analytic  solutions  do  exist  for  cer¬ 
tain  special  cases,  in  particular  the  one  dimensional 
case  which  we  consider  below. 

The  one  dimensional  form  of  eqn.  (6.1)  w ith  the 
constant  Y  set  equal  to  zero  is 


d  _  rvA  ic  l  <_p  _  Cj.  y  (_p  3 


(6.2) 


<« 


' in  ‘  ’ r  . 


0  f 


o  , 


rn 


■) 


o  r 


. 


n'Z::w 


' 


.  < 


‘"v  '  ‘  0  '  :  •'  ‘ 


■  • 


0 


. 


.  ■ 


d'mjJ'gnos 


36. 


The  following  substitutions  are  made  into  the  equation: 

L 9  +  u  = 

T&  A 


A 


A 

C 


(6.3) 


The  result  is 


(6.4) 


where  ^  and  "u”  are  now  dimensionless. 

Integrating  with  respect  to  l|^  *  we  obtain  the 

first  integral  of  eqn.  (6.4),  namely, 


(6.5) 


This  equation  can  be  solved  in  terms  of  elliptic  functions. 
The  form  of  the  solution  depends  on  the  constant  of  inte¬ 
gration  C.  The  solutions  of  this  equation  are 

t  i )  C  ~  o  ,  4>  - 

cos  h  (.  u  +  u-o) 

(J)  —  _2_  1 1  <c~ _ 

A6  COsU( 


(6.6) 
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K1  C  r°i  +  c  -/•  3 
-3.  f  <\  +c. 


CP  ^  J_  j?  ^  dn 

{Tu^i 


(  X  4  x  o  ^ 

v  A  jThJ^r  ) 


(6.7) 


o  >0-9  ^  ^  7  =  12  *c 

3  4  yr+c 


X  4  X  o 

A  /  z-  fc* 


(6.8) 


CivJ  C  -  9  ,  No  real  solution. 

Here  "k"  is  the  modulus  of  the  elliptic  function.  The 
constants  "k"  and  "xoM  are  arbitrary  constants  of  inte¬ 
gration. 

All  of  these  solutions  are  finite  everywhere.  The 
first  solution  has  clearly  decaying  exponential  asymptetic 
behaviour.  The  second  and  third  solutions  are  periodic  func¬ 
tions  whose  period  and  amplitude  depend  on  the  modulus  "k". 

The  three  dimensional  case  with  spherical  symmetry 
is  now  considered.  Eqn.  (6.1)  then  becomes 
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eqn.  (6.9)  becomes 


Now  the  substitution, 


(6.11) 


(6.12) 


is  made, 


where  "A"  is  a  constant.  The  result  is 


The  constant  nA"  is  chosen  so  that  the  relation, 


A-  A.5  + 

is  satisfied.  The  resulting  equation  is 


(6.13) 
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(6.14) 

If  the  function  "f"  goes  to  zero  sufficiently 
rapidly  at  large  distances  from  the  origin,  the  linear 
term  in  "f"  on  the  right  hand  side  of  the  equation  becomes 
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the  dominant  term  in  the  differential  equation.  Therefore 
if  the  constant  "A"  is  less  than  one  in  absolute  value, 
the  function  "f"  goes  as  exp  (  -  Jl  -  A2  x) . 

Only  two  cases  of  eqn. (6.14)  have  been  investigated. 
These  are  A  - ± 1  and  A  -  0.  The  second  case  will  be  dis¬ 
cussed  more  fully  later.  For  the  first  case,  eqn.  (6.13) 
gives 


_  2. 
+  3 


Eqn.  (6.14)  then  becomes 
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If  the  function  Mf"  is  written  as  a  series  of  the 
inverse  powers  of  nxM  and  substituted  into  the  equation, 
an  algebraic  solution  can  be  obtained  by  equating  powers 
of  "x**1" .  This  solution  is 


I  3  x 


which  is  finite  at  the  origin  and  goes  as  the  inverse 
power  of  Mx"  asymptotically.  The  function,  ^  is  thus 
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which  is  finite  and  non-zero  at  the  origin  and  at 
infinity. 

Solutions  exist  for  A2<  1  and  A  1  =  1,  which 

have  characteristically  different  asymptotic  behaviours; 
exponential  in  the  former  and  algebraic  in  the  latter. 
This  means  that  the  value  of  the  potential  at  infinity, 
"A",  greatly  influences  the  nature  of  the  solutions. 
This  clearly  illustrates  the  importance  of  the  choice 
of  gauge  in  the  original  time  dependent  equation. 

If  A  =  o  in  eqn.  (6.14),  the  equation  becomes 


(6.18) 


If  the  function,  Lp  ,  is  to  be  finite  everywhere ,  the 
function,  f,  must  be  zero  at  the  origin  and  finite  at 
infinity.  Since  at  large  distances  from  the  origin  the 
first  term  on  the  right  hand  side  is  the  dominant  one, 
the  asymptotic  form  of  "f"  is  exp(-x) .  If  the  substi¬ 
tution  f =  -f  is  made  into  the  equation,  it  is  unchanged. 
Therefore  if  "f"  is  a  solution  to  the  equation,  "-f"  is 
also  a  solution. 


If  a  change  of  variables  is  made  to 
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the  equation  is  transformed  to 


d2-f  ^  _i_  dj  =  1 

J  —  i  d  2  2 


A  3 


d**  *  4* 


3  Z 


(6.20) 


The  boundary  conditions  are  now  that 


•f  =  o  3  2-  =•  ° 


-f  -  O  ,  2  =■  1 


(6.21) 

We  also  have  that  near  the  origin,  "f"  is  linear  in  "z". 

Eqn.  (6.20)  has  been  solved  numerically  on  the 
LGP-30  computer  at  the  University  of  Alberta.  The  solution 
was  carried  out  by  using  a  step  by  step  forward  integration. 
The  following  approximations  were  made  for  the  values  of 
the  derivatives  at  any  one  point: 


I  h  +  l  “  -P  n  -1 


(6.22) 


2  In 


Here  Mh"  is  the  length  of  the  interval  between  successive 
points.  These  formulas  are  substituted  into  eqn.  (6.20) 
and  the  resulting  equation  is  solved  for  f  , .  The 
result  is 


42. 


To  find  the  value  of  the  function  at  any  one  point 
by  this  method,  the  value  of  the  function  at  two  preceeding 
points  must  be  known.  The  function  was  set  dqual  to  zero 
at  the  origin  and  the  value  of  the  function  at  the  first 
point  was  set  arbitrarily.  Then  the  value  of  the  function 
was  calculated  at  each  point  with  the  use  of  the  above 
formula,  up  to  the  point  z  «  1.  Since  the  boundary  condition 
requires  that  the  function  be  zero  at  this  point,  the  initial 
value  was  varied  until  the  final  value  became  zero. 

By  this  method,  three  solutions  have  been  obtained 
for  eqn.  (6.20)  which  fit  the  boundary  conditions.  The 
solutions  are  characterized  by  the  number  of  nodes  which 
they  contain.  Thus  the  first  solution  -  f  /x  has  no 
nodes,  ,  has  one  node  and  has  two  nodes.  Pre¬ 

sumably  there  are  many  more  solutions;  perhaps  an  infinite 
number.  More  solutions  could  not  be  obtained  because  of  the 
limited  accuracy  of  the  machine. 

The  functions  have  been  calculated  with  16,  128,  and 
512  points  in  the  interval  between  0  and  1.  Tables  I  and 
II  give  "f"  as  a  function  of  "z"  as  calculated  on  the 
computer.  Only  32  points  are  listed  for  each  solution. 
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16-point 

128-point 

512-point 

z 

solution 

solution 

solution 

0.00000 

0.000000 

0.000000 

0.000000 

0.03125 

• 

0.146789 

0.146917 

0.06250 

0.287915 

0.293504 

0.293757 

0.09375 

0.439997 

0.440370 

0.12500 

0.575140 

0.586038 

0.586523 

0.15625 

0.731300 

0.731883 

0.18750 

0.859569 

0. 8753^3 

0.876012 

0.21875 

1.017598 

I.OI8338 

0.25000 

1.137402 

1.157341 

1.158129 

0.28125 

1.293676 

1.294488 

0.31250 

1.402504 

1.425507 

1.426316 

0.34375 

1.551521 

1.552297 

0.37500 

1.645690 

1.670168 

1.670884 

0.40625 

1.779652 

1.780275 

0.43750 

1.854061 

1.877928 

1.878421 

0.46875 

1.962723 

1.963051 

0.50000 

2.010752 

2.031570 

2.031695 

0.53125 

2.O8I874 

2.081761 

0.56250 

2.095632 

2.111011 

2.110628 

0.59375 

2.116463 

2.115789 

0.62599 

2.087702 

2.095976 

2.095005 

0.65625 

2.047747 

2.046493 

0.68750 

1.969645 

1.970621 

1.969127 

0.71875 

1.864271 

1.862607 

0.75000 

1.734007 

1.729344 

1.727586 

0.78125 

1.567532 

1.565764 

0.81250 

1.388784 

1.381567 

1.379875 

0.84375 

1.175114 

1.173579 

0.87500 

0.959046 

0.952576 

0.951270 

0.90625 

0.718831 

0.717815 

0.93750 

0.482435 

0.478931 

0.478246 

0.96875 

0.237806 

0.237473 

1.00000 

0.000000 

0.000000 

-0.000021 

TABLE  I 


FIRST  SOLUTION  OF  EQUATION  (6.20) 
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fl 

f2 

z 

128-point 

128-point 

0.00000 

0.000000 

0.000000 

0.03125 

0.922109 

4.316718 

0.06250 

1.826049 

7.001103 

0.09375 

2.677376 

7.211515 

0.12500 

3.429140 

5.971093 

0.18750 

4.030550 

4.398816 

0.18750 

4.439408 

2.954004 

0.21875 

4.634133 

1.721698 

0.25000 

4.619803 

0.670752 

0.28125 

4.425375 

-0.247380 

0.31250 

4.094164 

-1.071010 

0.34375 

3.672684 

-1.822504 

0.37500 

3.202162 

-2.507862 

0.40625 

2.714333 

-3.117395 

0.43750 

2.230828 

-3.627526 

0.46875 

1.764585 

-4.005H6 

0.50000 

1.321993 

-4.215533 

0.53125 

0.904953 

-4.233714 

0.56250 

0.512582 

-4.054591 

0.59375 

0.142488 

-3.697656 

0.62500 

-0.208254 

-3.202726 

0.65625 

-0.542297 

-2.618971 

0.68750 

0.71875 

-0.861185 

-1.164283 

-1.992823 

-1.395654 

0.75000 

-1.447212 

-0.740834 

0.78125 

-1.699488 

-0.145099 

0.81250 

-1.901250 

0.427225 

0.84375 

-2.019615 

0.976756 

0.87500 

-2.007210 

1.489378 

0.90625 

-1.809018 

1.897868 

0.93750 

-1.385759 

1.999677 

0.96875 

-0.752024 

1.407735 

1.00000 

— - 

-0.000001 

0.000006 

TABLE  II 


SECOND  AND  THIRD  SOLUTIONS  OF  EQUATION  (6.20) 
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Figures  I  -  VI  are  graphs  of  <-p  and  "f"  as  functions  of  "x"  . 

When  the  number  of  intervals  is  increased,  the 
value  of  the  function  tends  to  increase  near  the  origin 
and  decrease  near  the  point  z  =  1.  The  first  solution  is 
the  only  one  which  has  been  calculated  at  all  three  inter¬ 
vals.  The  maximum  percentage  difference  between  the  16- 
point  solution  and  the  128-point  solution  is  about  2$, 
while  that  between  the  128-point  solution  and  512-point 
solution  is  about  0.1$.  The  128-point  solutions  were  used 
in  all  the  calculations. 

Maxwell's  equation  for  the  static  field  is 

dw  E  =  4np 


where  p  is  the  charge  density.  The  right  hand  side  of 
eqn.  (6.1)  is  identified  with  the  charge  density  for 
the  particle 
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i£-3 
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(6.24) 


The  total  charge  of  the  particle,  Q,  is 


a-r 


<L  V  =  - 


4  71 


V  z  d  V 


Using  the  divergence  theorem  we  have 


J'  V  4?  -  d  If 

Since  goes  exponentially  at  large  distances 


Q  -  -  -L 

4T7 


from  the  origin,  the  surface  integral  over  a  sphere  goes 
to  zero  as  the  radius  of  the  sphere  goes  to  infinity. 
Therpfore  the  total  charge  of  the  particle  is  zero  and 
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this  equation  must  describe  a  neutral  particle.  Since 
we  are  dealing  with  a  static  case,  no  intrinsic  magnetic 
moment  or  spin  is  involved. 

The  potential  for  a  charged  particle  would  have 
the  form  lP  -  A  +  q/r  asymptotically,  where  "A"  is  a  con¬ 
stant  and  "q"  is  the  total  charge  of  the  particle.  No 
solution  of  eqn.  (6.1)  has  this  asymptotic  behaviour. 

Thus  a  spherically  symmetric  solution  of  this  equation  could 
not  describe  a  charged  particle.  However,  since  charged 
particles  are  described  quantum  mechanically  by  complex 
fields,  it  is  possible  that  a  description  of  a  charged 
particle  could  be  obtained  in  our  static  case  by  carrying 
out  a  gauge  transformation,  in  which  case  according  to 
eqn.  (5.12),  E  would  be  replaced  by  the  complex  function 
Eexp( igS/ftc) . 

Figures  VII  -  IK  are  plots  of  p  and  rlp  as  functions 
of  the  radial  distance,  r,  for  the  three  solutions.  The 
graphs  of  the  charge  density  for  the  higher  solutions 
exhibit  an  orbital  arrangement  with  bands  of  concentrated 
charge  separated  by  regions  with  little  charge.  If  more 
solutions  exist  for  this  differential  equation,  one  would 
expect  them  to  show  more  of  these  charge  bands.  The  first 
solution  has  tv/o  of  these  charge  bands  of  opposite  sign. 

The  total  charge  in  one  of  these  hands,  q,  is 


q  =  3.016  g/G 


(6.25) 
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The  accuracy  of  the  first  machine  solution  has  been 
checked  by  various  means.  The  first  criterion  is  that 
the  integral  of  the  charge  density  over  all  space  is 
zero . 


(6.26) 


The  integrals  are  given  in  Table  IV.  The  difference  in 
the  integrals  for  the  16-point  solution  is  about  0.09%, 
and  the  difference  for  the  128-point  solution  is  about 
0.006%.  This  seems  to  indicate  that  the  machine  solutions 
are  quite  accurate,  at  least  for  the  first  solution. 

Another  attempt  to  obtain  some  indication  of  the 
existence  and  accuracy  of  the  solutions  is  given  below. 

Equation  (6.11)  can  be  treated  as  an  inhomogenous 
linear  differential  equation  with  the  cubic  term  as  the 
source  function. 


P  -- 

3 


(6.27) 


The  Green’s  function  for  this  equation  is 


(6.28) 


'  ^ 


V 


•  JV 


1  " 


'■')  ■’  0 


'I  V  ’  'I  * 


;■  -ro 


"\ 


59. 


Therefore 
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Integration  with  respect  to  the  polar  angle  gives 
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fo* 


An  analytic  approximation  was  made  for  the  function 
This  approximation  was  of  the  form 

h  =  A  Ce-x-£'bxJ 

^  -  4.1  3  b  «  3  ■ 


(6.30) 


A  3 


If  this  is  substituted  into  eqn.  (6.29),  the  result  is 

(c  -ffi  r  (o  -  2  66  5-74  )  d"* 

+  e'^fz  El  (u)  -13.35EcC4v4  6-xJ  f2fl.7Ei  C^-9x) 

-  ^,3$  e  c  ( 4, 35-x;  J  (6.31) 

~  E  4-  FcC4x)  —l  ^  .  35"  Ec^fe^S-x)  +  2i.7Et( 

-  II.3.S-  Ec  Cll.?S^)J 


where 


E^U)=  Xc^t 


C?D  t 

el  dt 


(6.32) 
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The  functions  f^.  and  are  tabulated  in  Table  III. 

An  improved  approximation  was  made  of  the  form 

+-t  -  A  C  -  e  -  4.x  e  ) 


A  -  4 . 7  ,  t»  -  4  )  acl.4-  ,  C-7 


(6.33) 


The  result  after  substitution  into  eqn.  (6.29)  is 


A3 


(0.18032,0)  £  'x-(o.o  s'as'jc  "/r^4.  (0.0  3440^/  e*4  *  * 
4  (o.QOlMOl  -  a.  0/2473*  )  Q.’Zni 

+  -e._x  La  E: C  ax)  -  <5  E «:  (sx)  +  as.i  Et(a> J 

-11.4  EcCiucJ  4  4.2  EcCi4x)J 


*  *  L-4  E  u  C4x)  -  2  I  frC  C7X  )  +  3  4.  2  E  i  Cl  o*) 
-  2  1,4  Ft  Cl3X  )  -t-  4 .  2  EcClhxJJ 


(6.34) 


These  functions  are  tabulated  in  Table  III. 
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X 

fo 

computer 

FIRST 

APPROXIMATION 

SECOND 

APPROXIMATION 

fG 

ft 

fG 

0.00 

0.0000 

0.000 

0.000 

0.0000 

0.0000 

0.10 

0.7296 

0.924 

0.885 

0.7755 

0.6792 

0.22 

1.4434 

1.572 

1.630 

1.5120 

1.5052 

0.36 

1.9392 

1.920 

2.055 

1.9749 

1.9899 

0.52 

2.1163 

2.013 

2.156 

2.1173 

2.1439 

0.70 

2.0249 

1.914 

2.024 

2.0138 

2.0363 

0.92 

1.7535 

1.676 

1.744 

1.7448 

1.7570 

1.20 

1.3784 

1.341 

1.364 

1.3752 

1.3789 

1.60 

0.9414 

0.930 

0.930 

0.9410 

0.9409 

2.30 

0.4705 

0.470 

0.465 

0.4707 

0.4702 

TABLE  III 


TrST  OF  FIRST  SOLUTION  WITH  GREEN1 S  FUNCTION 
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It  was  found  that  the  maximum  deviation  between 
the  analytic  approximation,  (6.30)  ,  and  the  Green’s  function 
solution  is  about  7%  •  With  the  improved  approximation, 
(6.33)9  ^he  maximum  deviation  is  about  1.3 %•  These  results 
suggest  that  even  better  analytic  approximations  to  the 
machine  solution  would  probably  show  even  less  deviations 
from  the  Green's  function  solution. 

Since  a  Lagrangian  has  not  been  obtained  for  eqn. 
(5.13)9  lhe  energy-momentum  tensor  cannot  be  determined. 
Hence  the  total  energy  of  the  particle  cannot  be  evaluated. 
The  assumption  is  made  that  the  total  energy  of  the  system 
can  be  taken  as  the  space  integral  of  the  product  of  the 
charge  density  and  the  potential,  i.e. 


£  - 


p  d  V 


(6.35) 


where  £  is  the  total  energy.  From  eqn.  (6.24)  we  have 

£--■5*=  [(-%  -Jl1  JV 

471  I  V.  A2-  3  ' 
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(6.36) 


iWx 


These  integrals  have  been  calculated  for  the  three  solutions 


and  have  been  tabulated  in  Table  IV. 
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CHARGE  INTEGRALS  FOR  FIRST  SOLUTION 


n  oc 

J 

"^0 

Jo  3 

16-point 

solution 

4.29697 

4.30075 

128-point 

solution 

4.366975 

4.366726 

ENERGY  INTEGRALS 


jo°%Zx2^  X 

jo 

16-point 

first  solution 

4.422716 

53.67182 

128-point 

first  solution 

4.509968 

54.12047 

128-point 

second  solution 

28.36375 

340.6186 

128-point 

third  solution 

85.97097 

1031.239 

TABLE  IV 
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It  is  interesting  to  note  that  1^  =  12I2  for  each 
of  the  three  solutions  obtained.  Thus  if  the  energy  were 
some  linear  combination  of  the  two  integrals,  different 
from  (6.36),  the  ratio  of  the  energies  for  the  three 
particles  would  be  unchanged. 

For  the  first  solution  we  have 
I  z  •=•  4..  S' o  &  8 

■=  1  a .  o  4  0 1  (c 

'b 

This  gives  a  value  for  the  energy  of  the  first  particle, 


£0  ^  /  3  .  S  3  me  *- 

G 


(6.37) 


If  the  mass  "m"  in  the  equations  is  to  be  taken 

as  that  of  the  lowest  state,  i.e.  the  first  particle, 

the  energy  from  the  first  solution  must  be  identically 

p 

equal  to  me  .  This  gives  a  value  of  13.53  for  the  coupling 
constant,  g2/fic,  which  is  of  the  order  of  the  coupling 
constant  for  the  "strong"  interactions. 

Since  the  value  of  the  coupling  constant  for  electro¬ 
magnetic  interactions,  o*  ,  is  known,  the  ratio  of  "g", 
the  basic  unit  of  charge  in  this  system,  to  "e"  the  charge 
on  the  electron  can  be  calculated.  We  get 
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and 
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e 

It  is  interesting  to  note  that  g^/e^  is  within  one  per 

cent  of  the  nucleon  electron  mass  ratio. 

We  had  earlier  in  eqn.  (6.25)  that  the  charge  in 
one  of  the  bands  in  the  first  solution,  q,  is  g/G. 

This  charge  is  then  equal  to 


q  -  0.2229  g  -  9.60  e 


(6.39) 


The  energies  of  the  second  and  third  particles  with 
respect  to  the  first  work  out  to  be 
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(b  .  a  9  3  2  0 

(6.40) 

1 9 .  0  &  ° 

(6.41) 

Interaction  potentials  between  the  particles  have 
been  calculated  using  analytic  approximations  to  the 
machine  solutions.  We  take  the  interaction  potential 
V-jj,  between  two  particles  as 
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Using  eqn.  (6.24)  and  the  divergence  theorem,  we  obtain 
the  following  result: 


j" Lpj  df' 


Therefore 


Since  for  a  given  solution  (-Pi  ,  -  <-Pi  is  also  a  solution, 

it  follows  that  the  interaction  potential  can  be  ±  V.  .. 

In  calculating  the  interaction  potential,  we  have 
assumed  that  the  charge  distribution  of  one  particle  is 
unchanged  in  the  field  of  another  particle.  Since  there 
will  be  mutual  distortion  when  the  particles  are  close 
together,  the  interaction  potential  is  only  qualitatively 
correct  near  the  origin.  In  doing  the  calculations,  it 
was  found  that  the  interaction  potential  could  not  be  deter¬ 
mined  accurately  near  the  origin.  There  would  be  a  further 
modification  if  the  two  particles  were  in  relative  motion 
since  the  interaction  is  velocity  dependent. 

The  interaction  potentials  have  been  calculated 
between  two  particles  of  the  first  kind  and  between  a 
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particle  of  the  first  kind  and  one  of  the  second  kind. 

This  was  done  by  making  analytic  approximations  for 

0  ^  t_p,  and  p<s  as  calculated  by  the  computer. 
Interaction  potentials  involving  the  third  particle  or 
the  one  between  two  particles  of  the  second  kind  could  not 
be  determined  accurately  because  of  the  difficulty  of 
approximating  the  potentials  and  charge  distributions 

sufficiently  closely. 

The  approximations  taken  were  of  the  form 
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These  potentials  are  tabulated  in  Table  V  and  plotted  in 
Figures  XII  and  XIII. 

From  eqn.  (6.43),  we  get  the  value  of  the  interac¬ 
tion  potential  at  the  origin  to  be 


If  the  two  particles  are  of  the  same  kind,  we  see  from 
the  definition  of  the  total  energy,  eqn.  (6.35)?  that  the 
interaction  potential  at  the  origin  is  equal  to  the 
total  energy  of  one  of  the  particles.  Thus  we  have  from 
eqn.  (6.37) 


Vo o Co)  =  I  3  .  & 3  r^c2- 

G 


(6.50) 


The  value  of  Vq^(o)  has  been  evaluated  by  numerical 
integration  of  the  functions  from  the  computer.  This  is 


Voi  M  =  1  ?  •  4  38  !21±Z 


G 


(6.5D 


We  see  that  these  values  are  quite  close  to  those  obtained 
using  the  analytic  approximations. 
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We  note  again  that  -Vqq  and  -V^  are  also  possible 
interaction  potentials,  since  the  sign  of  the  interaction 
potentials  depend  on  the  relative  signs  of  the  individual 
particles.  Thus,  for  example,  VQ0  can  be  essentially 
attractive  as  well  as  repulsive.  Therefore  V0q  in  Figure 
XII  is  shown  to  be  essentially  repulsive  because  the  signs 
for  the  potentials  of  the  interacting  particles  were  the 
same.  Whereas  if  they  had  been  taken  as  different,  the 
interaction  potential  would  have  be^n  essentially  attractive. 
The  same  arguments  apply  for  V^. 

If  we  choose  the  mass,  m,  to  be  the  mass  of  the 
neutral  pi-meson,  and  the  coupling  constant,  G,  equal 
to  13.53?  the  energy  scale  in  the  case  of  VQ0  would 

p 

correspond  to  a  value  at  the  origin  of  135  Mev,  i.e.  m-n-c  . 

In  the  case  of  Vq^?  the  minimum  would  correspond  to  an 
energy  of  64  Mev  or  0.47  m-n  c^.  In  both  cases  of  course, 
the  distance,  x,  would  be  in  the  units  of  the  pi-meson 
Compton  wavelength,  1.4  x  10'  '-3cm. 

The  shape  of  Vq-^  is  rather  interesting  because  it 
indicates  the  possibility  of  forming  a  bound  state  of 
particles  one  and  two.  The  binding  energy  would  lie 
between  0  and  64  Mev,  so  that  the  total  energy  of  the 
bound  system  including  the  energy  of  the  individual 
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O 

o 

> 

V01 

•  0.0 

13.38 

12.36 

0.1 

11.54 

0.2 

12.39 

9.36 

0.3 

6.47 

0.4 

10.14 

3.40 

0.5 

0.54 

0.6 

7.63 

-1.88 

0.7 

-3.76 

0.8 

5 .46 

-5.11 

0.9 

-5.97 

1.0 

3.74 

-6.42 

1.2 

2.47 

-6.40 

1.4 

-5.64 

1.5 

1.24 

1.6 

-4.60 

1.8 

-3.52 

2.0 

0.284 

-2.55 

2.2 

-1.75 

2.5 

-  0.028 

-0.875 

3.0 

-  0.103 

-0.089 

3.5 

-  0.101 

-0.191 

4.0 

-  0.079 

0.243 

4.5 

-  0.056 

G.212 

5.0 

0.162 

5.5 

0.114 

6 .0 

0.078 

TABLE  V 
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particles  would  like  in  the  range  6.p  to  7.3  ra  ti  c  . 
that  in  our  static  case  the  mass  ratio  of  the  compound 
particle  to  that  of  the  first  particle  is  approximately 
the  same  as  the  ratio  of  nucleon  mass  to  pi-meson  mass. 
Needless  to  say,  the  spin  of  our  compound  particle  is 
certainly  not 
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VII.  CONCLUSIONS 

We  have  seen  on  the  basis  of  a  rather  simple  and 
perhaps  naive  static  model,  that  it  was  possible  to  obtain 
a  number  of  discrete  states,  which  one  may  hope  would 
represent  in  a  more  complete  formulation,  a  spectrum  of 
particles.  Such  a  formulation  would  require  an  unambiguous 
definition  for  the  energy-momentum  tensor,  which  we  feel 
could  only  stem  from  a  first  order  boson  equation. 

With  respect  to  the  energy-momentum  tensor,  we  would 
like  to  offer  the  following  tentative  suggestionto  explain 
the  fact  &s  has  already  been  mentioned  on  page  53?  that  for 
the  three  discrete  solutions  obtained,  it  was  found  that 
j'1  c_p  4  d  V  -  izf'  ^  z  d  V  .  Nov;  we  have  not  been  able  to  show 
that  such  a  relation  is  an  automatic  consequence  of  the 
static  equation  for  ,  which  leads  us  to  believe  that  such 
a  result  may  be  of  more  fundamental  significance.  Consider, 
for  example,  a  static  Lagrangian  of  the  form 

cC=  <*•  [i  ( +  ^  g>1-  7k  *4J 

where  "a"  is  an  arbitrary  constant.  Disregarding  as  we 
must  the  relativistic  invariance  of  the  Lagrangian,  the 
above  Lagrangian  nevertheless  does  yield  the  static 
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equation  for  Q  .  It  is  easy  to  show  that 

-  [Yx*dv  -  [r,3Civ=  ft 

which  in  our  case  is  zero,  and  this  is  certainly  a 
necessary  requirement  for  a  stationary  particle.  (T^  of 
course  depends  on  the  choice  of  "a").  In  the  time  dependent 
case  the  symmetrization  of  the  energy-momentum  tensor  may 
lead  in  the  static  limit  to  our  expression  for  the  energy, 
i.e.  j'  cp  p  d  V  ?  and  to  for  ^T^dV  (®i=l,2,3)  • 

Finally  we  note  that  in  the  case  of  our  algebraic  solution, 
the  above  considerations  do  not  apply.  Hence  the  requirement 
for  particle-like  solutions,  such  as  fT**  av  =o 
may  distinguish  the  physically  acceptable  solutions. 
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-  electromagnetic  field  tensor 

-  magnetic  field  vector 

-  electric  field  vector 

-  vector  potential 

-  scalar  potential 

-  current  density 

-  charge  density 
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